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We present here the performance of a simultaneous dual-species matter-wave accelerometer for 
measuring the differential acceleration between two different atomic species (®^Rb and ®®Rb). We 
study the expression and the extraction of the differential phase from the interferometer output. 

The differential accelerometer reaches a short-term sensitivity of 1.23 x gjy/Hz limited by the 

detection noise and a resolution of 2 x 10“®(; after 11000 s, the highest reported thus far with a 
dual-species atom interferometer to our knowledge. Thanks to the simultaneous measurement, such 
resolution levels can still be achieved even with vibration levels up to 3 x corresponding to 

a common-mode vibration noise rejection ratio of 94 dB (rejection factor of 50 000). These results 
prove the ability of such atom sensors for realizing a quantum based test of the weak equivalence 
principle (WEP) at a level of r; ~ 10“® even with high vibration levels and a compact sensor. 


I. INTRODUCTION 

Light pulse atom interferometry [Il,i has allowed to 
probe, in an extremely sensitive and accurate way, the 
incidence of inertial forces on the motion of cold atoms. 
In recent decades, the very high performance of such sen¬ 
sors have been demonstrated for measuring t he g ravity 
acceleration MU, Earth’s gravity gradient [8|-[l^. rota¬ 
tions [HHil and gravity-field curvature (l^ . That is 
why they have been greatly dev elop ed for many applica¬ 
tions such as inertial navigation [l^, geophysics, geodesy, 
sub-surface exploration, and metrology (^. 

Many experiments and projects worldwide, also prove 
that quantum sensors appear to be very promising tools 
for exploring many aspects of fundamental physics as the 
determination of the fine-structure constant a [l^ and 
the Newtonian gravity constant G [Hill , the detec¬ 
tion of gravitational waves [l^ , the exploration of short- 
range forces [13,1m quantum based tests of the Weak 
Equivalence Principle (WEP) [23 - l23| . 

The WEP, also called the Universality of Free Fall 
(UFF), is one of the pillar of the Einstein’s theory of 
general relativity. This postulate states that all bodies, 
regardless of their internal composition, are affected by 
gravity in a universal way, in particular all test parti¬ 
cles at the alike space-time point in a given gravitational 
field will undergo the same acceleration. Today, different 
theories predict its violation [ 23 , Ull in the purpose of 
unifying general relativity with non-gravitational funda¬ 
mental interactions as described by the standard model 
in a quantum approach. Matter-wave tests thus offer new 
alternatives by their intrinsically quantum nature, radi¬ 
cally different from their classical counterparts [IJl , that 
can bring new constraints and bounds on the unifying 
theories. 

Presently, the state-of-the-art WEP test by atom in¬ 
terferometry is at a level of few 10“® [131 • hi order to 
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improve the sensitivity of atom accelerometers, several 
projects under development aim to compare the free fall 
of two different atomic species during few seconds, as 
the sensitivity scales quadratically with the interrogation 
time, in 10-m-tall atomic fountains [29l-l^. drop towers 


32|, sounding rockets, parabolic flights [33j and in space 
3J-|^|. These projects are developed in parallel with 
works on the increase of the momentum transfer l37H40ll 


to enhance the enclosed area and further improve the in¬ 
terferometers sensitivities. The use of Bose Einstein con¬ 
densates as atomic sources is also prospected [H, |4l|, H^l . 

Toda y, t he sensitivity of most state-of-the-art gravime¬ 
ters [j, [dJl is limited by vibration noise. To avoid this 
effect in quantum based test of the WEP, a special at¬ 
tention must be paid to the simultaneous interrogation 
of both quantum proof bodies in order to benefit from 
an efficient common-mode noise rejection allowed by the 
differential measurement [23j |. 

In this paper, we present the characterization of a si¬ 
multaneous dual-species atom interferometer which si¬ 
multaneously handles both isotopes of rubidium (®^Rb 
and ®®Rb). This sensor allows us to study the speci¬ 
ficities of a differential acceleration measurement, and 
especially the rejection of the vibration noise which is 
currently at 94 dB. We also show that the achieved sen¬ 
sitivity of our experiment is compatible with a test of 
the WEP at a level of 2 x 10“®. This work is part of 
the context of quantum based tests of the WEP and con¬ 
tributes to the elaboration of future experiments, both 
ground and space, aiming to detect a WEP violation. 

The paper is organized as following: section [H] gives 
the formalism for calculating a complete expression of 
the differential phase in our dual-species Mach-Zehnder 
type interferometer. Section IIIII explains then how this 
differential phase is extracted from the elliptic interfero¬ 
metric signal. Section ITVl describes the experimental ap¬ 
paratus and the overall measurement sequence. Section 
IVl presents the improvements brought to the experiment 
leading to the new results on the differential accelera¬ 
tion measurement given in section IVII Finally, section 
En] deals about the common-mode vibration noise rejec- 
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tion and presents its theoretical and experimental limits. 
A brief discussion on a quantum test of the WEP is given 
as a conclusion in section IVIIIl 

II. DIFFERENTIAL PHASE EXPRESSION IN A 
DUAL-SPECIES ATOM INTERFEROMETER 

For simultaneously measuring the acceleration under¬ 
gone by ®^Rb and ®®Rb, a Mach-Zehnder type atom in¬ 
terferometer is realized. This interferometer consists of a 
sequence of three equally spaced light pulses, of duration 
T —2r —T, driving stimulated Raman transitions between 
the two fundamental hyperfine states of the atoms [dj. 
At the end, considering the isotope i, the proportion of 
atoms in each state depends sinusoidally on the phase 
difference between both paths of the interferometer A$i, 
being proportional to the acceleration Ui undergone 
by the isotope i along the Raman laser direction of prop¬ 
agation. To first approximation, by considering infinitely 
short Raman pulses {i.e. t T), the phase is given by 

i 

A*, ~ ■ s, - 2TO)r^ (1) 

where klg is the effective wave vector of the Raman laser 
associated to isotope i, a is the microwave chirp applied 
on the Raman frequency to compensate for the Doppler 
shift induced by gravity for both isotopes, and T is the 
time between two pulses of light. 

The two signals from the dual-species atom interferom¬ 
eter are sinusoidal functions of the interferometric phase 
and can be expressed as 

Ps7 = Psr - ^ cos(A$87) 

Ps5 = P85 - ^ COs(A$ 87 + (/>d) 

with 

(j)d = A$85 — A<i)87 (3) 


to each isotopes which are basically equal to the primi¬ 
tive integrals of the previous sensitivity functions. The 
interferometric phases can thus be expressed by following 
this formalism, 

A$87 = / f87■ 087(t) - 27ra) dt + 

A$85 = / f85(t)(k^g ■ 085 (t) “ 27ra) dt + 

with 087(t) = 0-1- d{t) and 085(t) = o -I- Ao -I- d{t). 
087,85 (i) are the accelerations of each isotope projected 
along k^g^^ relative to the inertial reference correspond¬ 
ing to the Raman mirror in our experiment. They are 
composed of terms of different origins: the mutual con¬ 
stant acceleration o (corresponding to the gravity accel¬ 
eration g for a ground gravimeter and to 0 for a whole 
instrument in perfect free fall), the additional accelera¬ 
tion of the Raman mirror d(t) and the potential WEP 
violation signal Aa (as dehned in Eq. (fTSll l. For its part, 
represent the phase shifts due to systematic effects. 

The calculation and the demonstration of the sensi¬ 
tivity function have already been carefully described in 
[4^, we have slightly extended the experimental frame¬ 
work wherein this function is calculated. We assume here 
a case where the Rabi pulsations fli, D2, D3, associated 
to each of the three pulses of the interferometer, are not 
equal one to another (Di ^ U ,2 ^ ^ 5 ) and do not realize 
a perfect 7r/2-7r-7r/2 sequence (Di^2,3'r ^ 7r/2). Experi¬ 
mentally, these differences in Rabi frequencies may come 
from many sources. For instance, they can be induced 
by Raman laser power fluctuations between each light 
pulse. The transverse motion of an atom may also be a 
source of Rabi pulsation differences because of the spa¬ 
tial inhomogeneity of the Raman laser beam intensity. In 
our experiment, these differences are mainly due to addi¬ 
tional laser lines which generate a spatial dependency of 
the Rabi frequency. This point will be discussed more in 
details in section El and Appendix [BJ In this framework, 
the expression of the sensitivity function of the interfer¬ 
ometer can be derived from the following definition : 


where ^87,85 are the proportions of atoms in the upper 
hyperfine ground state, ^87.85 ^he offsets of the pop¬ 
ulation measurements, (^87,85 are the fringe amplitudes, 
A$87 is the interferometric phase for ®^Rb as expressed 
in Eq. (HD and (j)d represents the differential phase between 
the two species. We are interested in measuring this last 
quantity which is directly linked to the differential accel¬ 
eration undergone by atoms. 

The sensitivity function formalism has been used 
in this paper in order to express the more strictly as pos¬ 
sible this differential phase. This formalism allows to take 
into account the finite duration of the Raman pulses. The 
sensitivity function gg (t) is linked to the step response of 
the atom interferometer to an infinitesimal variation Sip 
of the Raman laser phase. In this paper, we are inter¬ 
ested in the impulse response of the Mach-Zehnder type 
atom interferometer to accelerations. These responses 
are given by the response functions f 87 , 85 {t) associated 


glit) = 


Ci sin (A^i) 6ip^o Sip 


6P^ {dip,t) _ A$i((5:^,t) 


S(f^O 


Sip 


(5) 


From the sensitivity function, the interferometric 
phase Ad^i can be evaluated for arbitrary evolution of 
the phase of the Raman laser ip{t) 


= J glit) dipit) 



dt 


( 6 ) 


In the frame associated to the free falling atom, the 
phase of the Raman laser experienced by the atom is 
given by ip{t) = where r) (t) is the position of the 

atoms i compare to the inertial reference. The equation 
thus shows that the sensitivity function corresponds 
to the impulse response of the interferometer relative to 
the velocity of the free falling atom. 
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FIG. 1. Mach-Zehnder type Atom interferometer time di¬ 
agram and associated acceleration response functions f{t) 
and velocity response functions gs{t). The black dashed line 
curves correspond to fli = fl 2 = fls = ^, the red solid 
line curves correspond to flir = 0.851-, fl22r = O.OOtt and 
flar = 1.08|- similar to values of the Rabi frequencies associ¬ 
ated to ®®Rb in our experiment. Here, T = 3 ms and r = 1 
ms, the small value of T and high value of r are arbitrarily 
chosen for a good readability of the curve but have no con¬ 
nection with the reality of the experiment. 

The expression of the response function to acceleration 
is then given by: 

Mt) = - (7) 

The Fourrier transform of this function gives the accel¬ 
eration transfer function of the atom interferometer (see 
figure [71) . The calculation and the exact expressions of 
both gs{t) and functions are given in Appendix El 
Two response functions, in acceleration and velocity, of 
a Mach-Zehnder type atom interferometer, correspond¬ 


ing to two different sets of Rabi frequencies (fli = 172 = 
ria and 17i 7 ^ fl 2 Ha), are shown in FIGjTJ In the case 
where the Rabi frequencies perfectly match each other, 
the acceleration response function is rigorously even and 
is zero outside the light pulses sequence. This proper¬ 
ties illustrate the fact that the interferometer is sensitive 
to acceleration during the interferometric phase only. In 
the case where the Rabi frequencies are different one to 
another, the acceleration response function still remains 
even on the interval [—T — r, T -|- r] whereas the parity 
is broken between the first ([—T — 2t, —T — r]) and the 
third ([T -|- r, T + 2t]) pulses if Ha 7 Z Hi. It then follows 
that the function does not go to zero before the first Ra¬ 
man pulse, which is the signature of an interferometer 
sensitive to the atom initial velocity. When Ha > Hi, 
the acceleration response function is slightly above the 
response function corresponding to Hi = H 2 = Ha and is 
thus positive for t < —T — 2t {cf. FIG. [1]). 

Thanks to this formalism, the scale factor S impacting 
the measurement of a constant acceleration can be de¬ 
rived from Eq.([B]). For deriving S, we consider that the 
velocity of the atoms evolves as at. S is thus given by: 

■S' = 

— + 2T)('r + fl;-tan(^)-I-tan(^^)) 

( 8 ) 

This expression gives the measurement scale factor at 
first order in r and 1/H. An experimental validation of 
this expression will be presented in section IVII Gl and the 
complete expression of the scale factor is given in ap¬ 
pendix |^(c/. Eq. (IA4I) 1. It is important to notice that 
at first order, the scale factor depends on Rabi pulsa¬ 
tions of the first and third pulses only. It can be easily 
understood because of the interferometer symmetry. The 
atomic mirror efficiency will impact the contrast of the 
interferometer only whereas the atomic beam splitters ef¬ 
ficiencies will impact the contrast and the phase of the 
interferometer. 

We have now all the required tools to rigorously express 
the differential phase (pd- All the terms composing it are 
listed in Table [H They can be gathered in five categories. 

Eirstly, terms coming from a potential WEP violation 
signal. Term 1 is the WEP violation signal Aa impacted 
by the measurement scale factor. Terms 2 and 3 are 
the first order corrections respectively coming from the 
differences in wave vectors, Sk = — k^g, and response 

functions, Sf{t) = fssit) — fsrit), between both isotopes. 
Term 4 is the second order correction which combines 
both differences. 

Secondly, terms depending on the mutual constant ac¬ 
celeration a undergone by both species, which are not 
canceled because of potential differences between both 
scale factors. Term 5 corresponds to the differential phase 
shift induced by the constant acceleration because of the 
difference in wave-vectors and term 7 combines this ef¬ 
fect with the difference in response functions. A differ- 
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TABLE I. Phase terms composing the expression of the differential phase 4>d and gathered in five categories: WEP violation 
signal terms (Aa), mutual constant acceleration terms (a), vibration terms (a), velocity terms (w), and systematics. The second 
column presents phase shift terms according to the response and sensitivity functions formalism. The third column presents 
the dominant term of the expansion in the low frequency limit (cu ^ 0). The notations are mainly described in the text. For 
terms 1 to 10 the mismatch between fli and fls has been neglected and we consider that fli = fl 2 = fls- This mismatch 
can not be neglected for velocity dependent terms (11 to 14) because they only appear for fli 7 ^ fls, where refers to 

tan(^^) — jA tan(^^) The fourth column gives the numerical size of the phase shift terms according to our experimental 

conditions (T = 47 ms, r = 4 /rs, a = g = 9.81 m.s“^, Aa = 2 x m.s“^, a = 3 x 10“^ m.s“^, Av — 6 mrn.s”^, Vs 5 = 0.6 

m.s“^, 5k/kgff = 5 X 10“®, | f Sf{t) dt\/\ f fsiif) dt| = 3 x 10“®, 50/17 = 0.1 (between the first and the third Raman pulse and 
between both isotopes)). 


Term 

Phase Shift 

Dominant Term 

Size 



(DC term, w —>■ 0) 

(rad) 

A a Terms: 




1 

ktffAa f fs7{t) dt 

^effAa(r -|- 2 r)(r + j^tan(A^)) 

7 X 10"^^ 

2 

5k Aa f fs 7 {t)dt 

SkAaiT + 2r)(T + tan(fl^)) 

3.5 X 10"^® 

3 

klsAa j Sf it) dt 

^effAa(T -I- 2r)(j^ tan(Aj^) — tan(Aj^)) 

2 X 10"^® 

4 

SkAa f 5f(t) dt 

SkAaiT + 2 r)(jj|^ tan (^2 ^tan (^2 ^)) 

1 X 10"^° 

a Terms: 




5 

Ska f fs 7 {t)dt 

Ska(T + 2t)(T + tan( ^ 2 ^)) 

1.73 

6 

{klla - 2170) f Sf(t) dt 

(^eff® “ 27ra)(r -b 2 r)(j^ tan(Aj^) — tan (^2 ^)) 

1 X 10“’’ 

7 

Ska f 5 fit) dt 

SkaiT + 2 r)(jj| 5 - tan(Aj-^) — tan (^2 ^)) 

4.8 X 10”® 

d Terms: 




8 

5k f f87(t)dit) dt 

Skd(T + 2t)(T’ + tan( ^ 2 '^)) 

5.3 X 10“* 

9 

kis f 5fit)a(t) dt 

kl^aiT + 2r)i^ tan(fl^) - ^ tan(fl^)) 

3 X 10“® 

10 

5k f 5f{t)dit) dt 

SkaiT + 2 r)(jj| 5 - tan(fA-^) — tan (^2 ^)) 

1.6 X 10“’® 

V Terms 




11 

kesAv J gf (t) dt 

1.87 A.. sn87 

1.2 X 10“® 

12 

{kgffVss + 2ivatr) f Sgait) dt 

(keffVSS + ZTratr) ( (^87)!^ ) 

1.2 X 10“® 

13 

Skvss f 5 ^( 7 ) df 

Skvss 

6.1 X 10“® 

14 

5kv85 f Sgait) dt 

5kV8sij§^ - (^) 

6.1 X 10“’’ 

Systematics: 




15 

4>fE - 

4>se - 4>se 

- 


ence in response function coupled to a middling correc¬ 
tion of the Doppler shift by the Raman frequency chirp 
(i.e. 27ra 7 ^ induces an additional phase shift (term 

6 ). 

Thirdly, terms depending on vibration noise a(t) and 
which will cause some limitations of the common-mode 
vibration noise rejection. These limitations result from 
5k (term 8 ), 5f{t) (term 9) or both (term 10). 

Fourthly, when Di 7 ^ Da the interferometer becomes 
sensitive to a constant velocity of the atom along the Ra¬ 
man laser (reflected by the fact that gl{t) dt 7 ^ 0 ). 


This induces a differential phase shift if the vertical con¬ 
stant velocities associated to each isotope are not equal 
Av = Vs 5 — Vs 7 ^ 0 (term 11). It then exists some 
correcting terms coming from the difference in wave 
vectors (term 13 and 14) and in sensitivity functions 
5gs{t) = gs^{t) — gs'^it) (term 12 and 14). Term 12 de¬ 
pends also on the start time tr of the Raman chirp which 
is set as to correspond to the time when the atoms are 
released from the trap. 

Fifthly, it obviously exists some systematics (term 15) 
which must be taken into account and corrected. In our 
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experiment, the main systematics are the impact of addi¬ 
tional laser lines, the two-photon light shift, the Coriolis 
effect and wave-front aberrations of the Raman laser . 


From this equation, the common interferometric phase 
A$ 87 j can be “directly” eliminated by using inverse 
functions (arcsin and arccos). By using some algebra, 
each couple {xj,yj) leads to two solutions of (pd equal to 


III. DIFFERENTIAL PHASE EXTRACTION 

The two signals from the dual-species atom interfer¬ 
ometer are sinusoidal functions and thus parametrically 
describe an ellipse (c/. Eq. ([2])). In order to benefit from 
the coupling between both sensors few methods for de¬ 
riving the differential phase from this ellipse have already 
been developed. The operating range of a simultaneous 
dual-species atom interferometer is much larger than for a 
“standard” single-species atom gravimeter. Indeed, even 
when the acceleration fluctuations are greater than one 
fringe spacing, the ellipse remains visible and the differ¬ 
ential acceleration can still be derived from it. This is the 
signature of the correlation between both interferometric 
signals from both isotopes which is made possible by the 
simultaneous aspect of the measurement. 

A first solution, historically developed for gradiome- 
ters, consists in executing a least-square fitting method 
. This method permits a rapid extraction of the differ¬ 
ential phase but is not bias-free in the presence of noises. 
A Bayesian analysis is a much more comprehensive ap¬ 
proach for estimating the differential phase [U and 
will lead to an optimal estimator with negligible system¬ 
atic error. This method can also be employed in cases 
where the scale factor associated to each species are dif¬ 
ferent [ 4 ^. Nevertheless, this method requires an accu¬ 
rate statistical model for the interferometer parameters 
and is computationally intensive to implement. For our 
data processing we have used an alternative method re¬ 
ferred as “Direct Phase Extraction” that is described 
in this section. 

In our case, the ellipse parameter corresponds to the 
common interferometric phase A<I >87 shared by both iso¬ 
topes which can be swept by scanning the Raman fre¬ 
quency ramp a or by introducing some vibration noise. 
The differential phase (pd is related to the ellipticity, when 
(pd = 7''/2 the normalized ellipse is a circle, when (pd = 0 
the ellipse collapses to a line. More specifically, these 
curves are not closed Lissajous curves. Indeed, the fre¬ 
quencies of the two sine functions are slightly different be¬ 
cause of the different scale factors related to each atomic 
species. Nevertheless this aspect is clearly negligible in 
our case, considering the slight difference in scale factors 
and the amplitude of vibrations. The parametric curves 
are thus approximated by ellipses in the paper. After 
each interferometric cycle, a couple of points {xj,yj) is 
obtained. The measurement is re-iterated N times and 
N couples are finally obtained to describe the ellipse, such 
as: 


r x--2^ 

\ Xj — ^ I 

I _ o 
[ 2/i — 2 


= cos(A$ 87 y) 

= C0s(A$87j + (pd,j) 


(9) 


VjG[l,iV] 

(pd,j = arccos (^X 0 j ± ^(1 - xj){l - yj)^ , G [0, tt] 

. 

The final value of (pd is extracted by using a maximum 
likelihood estimation over all (pd,j- Firstly, the phase esti¬ 
mation likelihood function, containing the 2N solutions, 
is fitted by a gaussian function to find its maximum in 
order to have a first estimation of pd- According to this 
first estimation, for all j, only the closest value pdj is 
then kept among the two solutions {cf. Ea. (fT0)l l. Indeed 
for each couple of point j, only one solution pdj gives 
the correct value of pd on the interval [0,7r], the other 
being randomly distributed on [0, tt] when the ellipse is 
randomly parametrized. Finally, pd is given by the mean 
value of the second phase estimation likelihood function, 
now containing only N values pd,j associated to the N 
measurements {xj,yj). The one-sigma resolution is given 
by the standard deviation of the likelihood function over 
'/N- In our experimental conditions the dominant noises 
are gaussian noises. The resulting likelihood function is 
thus symmetric and its mean value corresponds to the 
more likely value of the differential phase. 

Typical simultaneous dual-species fringes are reported 
in Fig. (Ha). Figure [UJb) shows the experimental ellipse 
obtained by plotting the interferometric signal from ®®Rb 
versus the one from ®^Rb. The phase estimation likeli¬ 
hood function {cf. Fig. [2l[c)) represents the differential 
phase distribution from this ellipse {pd ~ 1.7 rad). 

This differential phase extraction method is quick and 
easy to implement and chiefly allows for a bias-free dif¬ 
ferential phase extraction on a large range of pd- We 
performed computing simulations in order to verify the 
bias-free aspect of this method. For different value of 
pd G [0, tt] and different amplitudes of noises, pd is es¬ 
timated for a large number of noisy ellipses. This esti¬ 
mation method is bias-free on the range [Stt/IG, IStt/IG] 
for an offset noise ctp/P of 1%, an amplitude noise ac/C 
of 2% and a non common phase noise of = 30 mrad 
(all corresponding to noises three times larger than our 
experimental noises). This method requires an a priori 
knowledge of the fringe amplitudes {Ci) and population 
offsets (P°). These parameters can be easily deduced by 
fitting the sinusoidal fringes, scanned either by the Ra¬ 
man frequency ramp a in a low vibration environment 
or by the vibrations themselves, where in that case, the 
atom sensor is correlated with a mechanical accelerom¬ 
eter. The probability density function of the atom 
interferometer measurements can also give access to this 
parameters. A last method can be to simply derive the 
fringe amplitudes and offsets from the mean value and 
the standard deviation of the probability density function 
without any fitting method. No significant biases on the 
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FIG. 2. Simultaneous interferometric signals for ®^Rb and 
®®Rb, with T = 47 ms. From top to bottom, (a) interfer¬ 
ence fringes as a function of the microwave chirp a, in black 
squares from ®^Rb in red dots from ®®Rb. (b) Signal from 
®®Rb is plotted vs signal from ®^Rb drawing the ellipse fit¬ 
ted by the solid blue line, (c) is derived from the phase 
estimation likelihood function (red histogram). The fringes 
and the ellipse contain 2400 points, each point is acquired 
in 0.25 s corresponding to a global integration time of 600 s. 
The data were recorded with a single direction of the Raman 
wave-vector -|-/CefF. 


differential phase estimation appears for relative errors 
less than 5% on the fringe contrasts and 1% on the off¬ 
sets. All the conducted numerical simulation tests prove 
thus the high reliability and robustness of this phase ex¬ 
traction method. 

In a “standard” single-species atom gravimeter the 
phase resolution is proportional to the inverse of signal 
to noise ratio (SNR), i.e. it scales as SNR“^ (sj. With 
the present Direct Phase Extraction method, this depen¬ 
dency is statistically estimated at SNR“°-^. This be¬ 
havior is empirically estimated by performing numerical 
simulations: for typical fixed values of offset noise, am¬ 
plitude noise and non common phase noise, the SNR is 
modified by changing the fringes amplitudes. This phase 
extraction method is then less sensitive to acceleration 
fluctuations than the one used in a gravimeter [^. This 
is here explained by the fact that the signal is integrated 
over all possible interferometric phases [0, 27r] and not 
only in the neighborhood of n/l where the interferome¬ 
ter’s response is the most sensitive. 


IV. EXPERIMENTAL APPARATUS 

For simultaneously measuring the acceleration under¬ 
gone by ®^Rb and ®®Rb, both isotopes are first trapped 
and cooled thanks to the same laser beams and magnetic 
field gradients to form two spatially embedded magneto¬ 
optical traps (MOT). Both clouds are released from the 
trap and the accelerations are simultaneously measured 
by the same Mach-Zehnder type atom interferometer 
based on stimulated Raman transitions. 

The experimental setup is mainly derived from sm. 
The cold atoms are trapped and cooled at the top of 
a high vacuum chamber made of glass and are then 
dropped over a distance of ~ 4 cm limiting T to be 
smaller than 47 ms. The atoms are interrogated during 
their free fall by the Raman laser beam which is retro- 
reflected by a mirror representing the inertial reference 
for both isotopes. This mirror is attached to a mechan¬ 
ical accelerometer (nanometrics TITAN) to monitor its 
vibrations. The whole sensor head - containing the vac¬ 
uum chamber, the magnetic shield consisting of four lay¬ 
ers of mu-metal, the magnetic coils and the optics - is 
placed on a passive vibration isolation table (Minus-K). 
This table is itself mounted on an excitation table ac¬ 
tuated by piezoelectric transducers (PZT). The isolation 
table resonance is tuned at the excitation frequency in or¬ 
der to filter higher harmonics to obtain a clean excitation 
at a well defined frequency. 

The laser system for addressing both ®^Rb and ®®Rb 
is based on the frequency doubling of Telecom sources 
[5ll| . A distributed feedback (DFB) laser diode at 1560 
nm is amplified in a 5 W erbium-doped fiber amplifier 
(EDFA) and then frequency-doubled in a periodically 
poled lithium niobate (PPLN) crystal. The frequency of 
the laser is controlled thanks to a beat-note with a refer¬ 
ence laser locked on a rubidium transition. The needed 
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laser lines for the dual-species experiment are then syn¬ 
thesized thanks to a fibered electro-optic phase modula¬ 
tor at 1560 nm. Concerning the cooling and the detection 
stages, four laser lines are needed : cooling and repump¬ 
ing laser lines for both isotopes. In that cases, the carrier 
frequency is tuned on the ®^Rb cooling transition. The 
three other laser lines are generated by injecting three mi¬ 
crowave modulation frequencies (1.126 GHz, 2.915 GHz 
and 6.568 GHz) into the phase modulator as follows : 


— Jcarrier^ 

— fcarrier 1-126 GHz, 

— /carrier + 6.568 GHz, 

= fcarr^er + 1-126 GHz -t 2.915 GHz. 

( 11 ) 

The power of each line is controlled by adjusting the 
modulation depth of each microwave frequency. Dur¬ 
ing the cooling stage the experimental parameters are 
adjusted to obtain close trap features for both isotopes. 
Indeed both cooling line powers are set to be the same 
for both isotopes and the repumping power of ®^Rb is a 
little stronger than the one of ®^Rb, that compensates a 
higher depumping rate due to spontaneous emission be¬ 
cause of its tighter hyperfine structure. Approximately 
one-third of the global laser power is lost in additional 
modulation lines far from any atomic resonances. Phase 
modulation is also used for generating the laser lines dur¬ 
ing the interferometric sequence. Both Raman pairs are 
generated by directly injecting the Raman difference fre¬ 
quencies associated to each isotope (i.e., 6.834 GHz for 
®^Rb and 3.035 GHz for ®®Rb), making the carrier fre¬ 
quency common to both Raman pairs (cf. FIGl3]). This 
way of laser frequencies generation leads to the creation 
of additional laser lines that can induce destructive in¬ 
terferences of the transition probability by driving “par¬ 
asites” Raman transitions (cf. FIGlS]). The Raman 
pair corresponding to ®^Rb is red-detuned by 0.59 GHz 
with respect to the excited hyperfine state F' = 2 and 
therefore the one corresponding to ®®Rb is red-detuned 
by 1.86 GHz with respect to F' = 3. The power of these 
two pairs is adjusted to obtain Rabi frequencies of the 
two photons Raman transitions as identical as possible. 

With this setup, few 10® atoms are loaded from a back¬ 
ground vapor into MOTs in 250 ms. The atoms are then 
further cooled down in an optical molasses phase of 28 
ms leading to a temperature of 1.6 /iK for ®^Rb and 2.6 
/iK for ®®Rb. Additional t rap loss collisions due to inter¬ 
species atomic collisions do not exceed 10-15% in 
our case. These results confirm that the additional laser 
lines do not have any significant impact on the cooling 
efficiency. Atoms are furthermore selected in the Zee- 
man sub-level mp = 0 of the hyperfine ground state (i.e. 

= 1 for ®^Rb and F = 2 for ®®Rb) to remain insensi¬ 
tive to parasites magnetic fields. The selection is made 
thanks to a microwave 7r-pulse which drives the transi¬ 
tion between the two hyperfine ground states for atoms 


f87 

J cooling 
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J repumping 
J repumping 



<—►^^Rb Raman Pair < ►^^Rb parasitic Raman Pair 



4 




®^Rb ®®Rb 

FIG. 3. Schematic of the laser spectrum at 780 nm during 
the interferometric sequence with atomic level systems and 
Raman transitions for both isotopes. In solid lines the Ra¬ 
man pairs used for the atomic mirrors and beam splitters, in 
dashed lines the “parasite” Raman pairs due to phase mod¬ 
ulation generation, red color corresponding to ®®Rb and blue 
to «'^Rb. 


in the sub-level mp = 0 only. A slightly blue-detuned 
optical pulse removes any residual atomic population. 

During the free fall, the interferometric sequence oc¬ 
curs in a vertical uniform magnetic field of 28 mG. The 
sequence consists in three Raman laser pulses of dura¬ 
tions T — 2t — T ,with T = 4 /is, equally spaced in time 
by T = 47 ms and realizing the -Kj2-17-1:12 sequence of 
the Mach-Zehnder type interferometer (cf FIGH]). The 
Raman laser pulses couple at the same time the states 
\F = l,mp = 0 > to |F = 2,mp = 0 > for ®'^Rb and 
IF = 2,mp = 0 > to |F = 3,m_F = 0 > for ®®Rb. Ex¬ 
actly the same microwave chirp |a| ~ 25.143 MHz.s“^ is 










































applied to both Raman difference frequencies in order to 
compensate the time-dependent Doppler shift induced by 
gravity (c/. eq. ([l])). This chirp is synthesized with the 
same DDS (Digital Direct Synthesizer) for both isotopes. 

Finally, the atomic population repartition between the 
two coupled states is measured for each species by fluo¬ 
rescence detection. The atomic cloud is illuminated by 
a vertical resonant detection beam and the fluorescence 
is collected thanks to a collimation lenses system and a 
photodiode on the perpendicular direction. The atomic 
cloud is illuminated by two successive sequences of three 
light pulses of durations 1.5-0.05-1.5 ms. The first se¬ 
quence induces the fluorescence signal from ®^Rb atoms: 
the first pulse detects atoms in F = 2, the second one 
fully transfers atoms from F = 1 to F = 2, the third 
pulse is identical to the first one and detects atoms ini¬ 
tially in F = 1. During the first and third pulse the laser 
frequency is slightly blue-detuned (-1-0.88 T, where T is 
the natural line-width of the transition) and no mod¬ 
ulation is applied to detect and push away the atoms. 
During the repumping pulse (second one) the 6.568 GHz 
modulation is injected into the phase modulator to gen¬ 
erate the repumping line with a power adjusted in order 
to suppress as much as possible the carrier frequency. 
The second sequence of three pulses induces the fluores¬ 
cence signal from ®®Rb and is conceptually identical to 
the first one: the first pulse detects atoms in F = 3, the 
second one fully transfers atoms from F = 2 to F = 3, 
the third pulse is identical to the first one and detect 
atoms initially in F = 2. In practice, during the first and 
third pulse the 1.126 GHz modulation is injected into the 
phase modulator to generate the cycling transition with 
a blue detuning of 0.7 T and with a power allowing to 
cancel the carrier frequency and thus minimize a cross 
signal from ®^Rb atoms. During the repumping pulse 
(second one) the 1.126 and 2.915 GHz modulations are 
injected into phase modulator to generate the repump¬ 
ing line. In that case, both cycling transitions can not 
be canceled. The fluorescence of the background vapor 
is finally recorded and removed for both isotopes. With 
this detection scheme, a detection noise on the transition 
probability of crp = 0.0027 for ®^Rb, and ap = 0.0035 
for ®^Rb, is obtained. 

The whole sequence is performed at a repetition rate 
of 4 Hz. 


V. PREVIOUS RESULTS & IMPROVEMENTS 

Before giving the last results of the dual-species atom 
interferometer, we recapitulate the previous results pre¬ 
sented in and the improvements that have followed. 

The resolution on the differential phase was about 
3 X 10“®g after an integration time of 15 min correspond¬ 
ing to a sensitivity of I x lQ~^g/\/Hz. These perfor¬ 
mances were limited by the detection noise of the exper¬ 
iment. The fringe amplitude was about 22 % for ®^Rb 
and 7 % for ®®Rb mostly because of spurious impact of 


additional laser lines generated by modulation and of our 
detection scheme. Thanks to the simultaneous aspect of 
the experiment, we highlighted a common-mode vibra¬ 
tion noise rejection between both isotopes higher than 
55 dB (rejection factor of 550) limited by the detection 
noise. Moreover, we performed a quantum based test of 
the WEP at a level of few 10“^ limited by uncertainties 
over systematics. This experiment was the first demon¬ 
stration of a simultaneous dual-species atom interferom¬ 
eter. 

The first step to improve the experiment was to reduce 
the impact of additional laser lines generated by phase 
modulation. The spurious impact on the fringe visibility 
is briefly presented in appendix |b1 and a more compre¬ 
hensive theoretical paper [5^ studies the additional lines 
impacts more generally. The main point to emphasize is 
the spatial dependency of the probability amplitude be¬ 
tween the two states coupled by the light. The atomic 
beam splitters and mirrors efficiency directly depends on 
this probability amplitude (given by the Rabi frequency 
H). As atoms are in free fall, the two-photon Rabi fre¬ 
quency depends on the time when the light pulse occurs 
and thus the transfer efficiency depends on T. T is set at 
47 ms to achieve the best sensitivity allowed by our ex¬ 
perimental set-up and the retro-reflecting Raman mirror 
position has been adjusted to obtain Rabi frequencies 
for each pulse, and for each isotope, as close as possi¬ 
ble. For instance in the current experiment the different 
Rabi frequencies are = Tr/{7.7fis), = 7r/(8/rs), 
= 7r/(7.9/rs), = 7r/(9.4/rs), D®® = 7r/(8.I/xs), 

Hi® = 7r/(7.4/rs), that is why the pulse duration r = 4 
fj,s has been chosen for realizing a Mach-Zehnder type 
atom interferometer t: 12 — tt — tt 12 . 

The non-common phase noise between both isotopes 
was reduced by optimizing the micro-wave source used 
for generating the Raman pair of ®®Rb. Mainly, a lower 
noise frequency reference at 10 MHz was implemented in 
the optimized micro-wave source. 

The detection has also been improved in order to in¬ 
crease the SNR. The interferometric signal (cf. Eql^) is 
given by the proportion of atoms P in the upper hyper- 
fine ground state at the output of the interferometer : 


where N 2 is the number of atoms in the upper state and 
Ni in the lower one. The detection aims to accurately 
estimate the ratio P by counting the atoms in each state. 
We name here and the detected signal asso¬ 
ciated to each state. Our detection scheme shows non- 
linearities when Ni is measured {i.e. is not propor¬ 

tional to A^i). The interferometric fringes are not per¬ 
fectly sinusoidal any more and it results a bias on the 
determination of the differential phase from the ellipse. 
These non-linearities have been taken into account as in 
the following. 

The detection sequence has been explained in section 
CYl During the first detection pulse, when is N 2 mea- 
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sured, a part of these atoms, espN 2 , are transferred to 
the lower state because of the depumping rate due to 
spontaneous emission (0 < Cgp < 1). Moreover, some 
atoms, e 2 N 2 , are not perfectly pushed away from the de¬ 
tection area because of the finite duration of the pulse 
(0 < £2 < 1). These effects counterbalance each other 
(when £2 increases, £sp decreases and inversely) as a func¬ 
tion of the detection laser power. When the laser in¬ 
tensity is low compare to the saturation intensity of the 
detection transition, a very small part of atoms are lost 
because of spontaneous emission whereas a huge part of 
atoms are not pushed away from the detection area, and 
inversely at high laser intensity. 

During the second detection pulse, when atoms in the 
lower state are re-pumped in the upper state, it is im¬ 
possible to completely cancel the blue detuned detection 
light. This residual pushes away a part of the atoms, 
(1 — £i)iVi, from the detection area (0 < £i < 1). 

Consequently, during the third detection pulse, oc 
£ifVi -I- £i(£sp -I- € 2)^2 atoms will be detected instead of 
Ni. This effects must be evaluated corrected. The atomic 
population repartition can then be reconstruct by 


P = 


Ndet 


js^det 


+ (1 “ Esp ~ € 2 )^ 2 ^ 


(13) 


Interferometric fringes shown in Fig. [21(a) are corrected 
from these non-linearities. We measured for ®^Rb, re¬ 
spectively for ®®Rb, £1 = 0.75 and £sp -b £2 = 0.12, re¬ 
spectively £1 = 0.8 and £sp -b £2 = 0.3. Micro-wave power 
fluctuations give rise to temporal fluctuations of these 
non-linearity coefficients. These temporal variations are 
taken into account by fitting the interferometric fringes 
and extracting the parameters £1 and £sp -b £2 at different 
moments of the experiment. 

The SNR is now approximately five times better than 
before [ 2 ^, the fringe amplitudes being currently about 
40% for ®^Rb and 35% for This disparity is ex¬ 

plained by the impact of the additional laser lines which 
still remains a bit stronger for ®^Rb and by the visibility 
loss induced by the higher temperature of ®®Rb and the 
velocity selection of stimulated Raman transitions [^ . 


VI. RESOLUTION, SENSITIVITY AND LONG 
TERM STABILITY 

The interferometric signals shown in Fig. |21 associated 
to a single direction of A:eff, result from a measurement 
acquisition of 600 s at a repetition rate of 4 Hz. By 
extracting the differential phase from the ellipse as ex¬ 
plained in section m 4>d is derived with a resolution of 
1.5 mrad corresponding to a resolution of 5.10“®g and a 
sensitivity of 1.23xl0“^(? at 1 s, with our differential ac¬ 
celerometer by assuming that we are limited by a white 
noise. 

In order to study the long term stability of the instru¬ 
ment, we alternatively record at 4 Hz the signal for the 



FIG. 4. Allan deviations on the differential phase (j>d ex¬ 
pressed in g unit 4>d/ {k^gT'^) with g ~ 9.81 m.s“^. In black 
and red dots the Allan deviation corresponding to the two 
reversed directions of fceff. The Allan deviation of the averag¬ 
ing over these two directions (blue dots) and its asymptotic 
behavior (blue line) is fitted by 295 ng/yr. The best sensitiv¬ 
ity achieved with our differential accelerometer (green dots, 
corresponding to data points from FIG. O and its asymptotic 
behavior (green line) is fitted by 123 ng/yT. The first data 
point at r = 50 s is given by the time required to obtain 
enough points to extract the phase from the ellipse (here 100 
points). 


direction +kes and then for the opposite one —keg. A 
value of the differential phase (/)d is thus obtain every 0.5 
s. The direction of k^g is changed drop by drop to elimi¬ 
nate some systematic effects (mainly the one photon light 
shift and the first order Zeeman shift) whose sign does not 
change with keg. We split the atom interferometer data 
points into groups of 100 consecutive points as a mini¬ 
mum of 100 points is required to derive a confident value 
of (pd with our phase extraction method. The interfero¬ 
metric phase A $87 is randomly distributed over several 
fringes thanks to vibrations as the vibration isolation sys¬ 
tem was turned off during the data acquisition. For each 
group composed of 100 points, a value and a standard 
error are derived for (pd- For each group of points, the es¬ 
timation of detection non-linearities parameters as well as 
the estimations of the fringe amplitudes and population 
offsets are realized by fitting the interferometric fringes. 
This fit is achievable thanks to the correlation between 
the atomic and the mechanical accelerometer that allows 
to reconstruct the interferometric fringes. 

Figure [4] shows the Allan deviation on the differential 
phase pd, expressed in g unit {pd/{klggT'^)), where g 
refers to the gravity acceleration {g ~ 9.81 m.s“^), asso¬ 
ciated with three possible set of points: -\-keg points only, 
—fceff points only or averaging over the two opposite di¬ 
rections. This averaging allows to reject the long term 
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drifts induced by the one-photon light shift fluctuations, 
indeed the Allan deviation does not contain the bump 
around 3 hours induced by these fluctuations. 

These Allan deviation behaviors prove also that the 
short term sensitivity is limited by white noises, namely 
the detection noise and the non-common phase noise be¬ 
tween both isotopes. In our experiment, the detection 
noise is induced by the laser frequency noise during the 
detection. This noise is estimated to limit the sensitivity 
at a level of 1.4xl0“^g at 1 s. The other noise source is 
the non-common phase noise induced by the phase noise 
of the two different micro-wave sources used for generat¬ 
ing both Raman pairs by phase modulation. By experi¬ 
mentally measuring the power spectral density of phase 
noise associated to each source, we estimated that the 
impact of this noise is lower than 0.65x at 1 s. We 
finally estimated a global white noise level limiting the 
sensitivity to 1.55x10“^^ at 1 s. This value is in good 
agreement with the sensitivities measured in FIG. SI 

The best experimental sensitivity (I.23xI0“^5 at I s) 
was achieved with data points from FIG. S] for a single 
direction of fceff. The sensitivity after the averaging over 
the opposite directions of fceff (2.95xl0“^g at 1 s, c/. 
blue points in FIGS]) comes from an different set of data. 
During this measurement session, the detection was not 
fully optimized degrading the sensitivity. 

The results shown in FIG. SI demonstrate that our dif¬ 
ferential accelerometer has a sufficient resolution, by in¬ 
tegrating the signal during only 5 hours, for testing the 
WEP at a level of 2x10“®. 


VII. VIBRATION NOISE REJECTION 
A. Experimental Results 

In our experiment, the differential acceleration mea¬ 
surement is performed with a simultaneous interroga¬ 
tion of both atomic species falling in a common refer¬ 
ence frame. The simultaneous aspect allows us to benefit 
from an efficient common-mode noise rejection which is 
not possible when atoms are alternatively handled. 

The sensitivity of most state-of-the-art gravimeter is 
limited by the vibration noise whereas the differential 
atom sensors intrinsically allows to efficiently reject it. 
The rejection of this noise, in an atom sensor, has al¬ 
ready been studied in a gravity gradiometer Q handling 
two cold atomic clouds of the same species separated in 
space. Our experiment has demonstrated for the first 
time the rejection of vibration noise with two embed¬ 
ded clouds of ®^Rb and ®®Rb simultaneously submitted 
to the same light pulse interferometric sequence. In our 
case, this kind of rejection can be extended to other en¬ 
vironmental perturbations, such as gravity gradients or 
rotations, which depend on the overlap of the two atomic 
ensembles. 

The vibration noise rejection is a critical point for a 
sensor which might be used on a moving platform or in a 



0.3 0.4 0.5 0.6 0.7 

Fraction of ®^Rb in upper state 


vibrations of 2.4x10'^g 



0.3 0.4 0.5 0.6 0.7 

Fraction of ®^Rb in upper state 


FIG. 5. Three ellipses of 2000 points plotted for three dif¬ 
ferent vibration amplitudes: 3.2 x 10~^g, 7.6 x W~^g and 
2.4 X 10~^g. The ellipses are progressively blurred because of 
additional inertial effects induced by the excitation (cf. sec¬ 
tion |VnDl). 


highly noisy environment. In the more specific context of 
a WEP test, it is necessary that the vibration noise does 
not limit the sensitivity of the instrument, i.e. terms 8, 
9, 10 presented in Table U must remain lower than the 
targeted one shot resolution. So the vibration noise re¬ 
jection ratio must be evaluated in order to determine the 
vibration level limit at which the instrument still benefit 
from an adequate sensitivity. 

We characterized the rejection of vibrations by shak- 
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Vibration ampiitude (g) 

FIG. 6. Resolution on the differential acceleration versus the 
amplitude of vibrations. The resolution is estimated from the 
ellipses, containing N — 2000 points, for different amplitudes 
of vibrations at 2.08 Hz. The correlation with a mechanical 
accelerometer allows to measure the vibrations, and to esti¬ 
mate the fringe amplitudes and the population offsets. These 
two parameters depends on inertial effects (cf. IVIIDll . these 
dependency is (red dots) or is not (black squares) taken into 
account. In dashed lines, the resolution is plotted according 
different values of the rejection ratio r = 88.2 dB (black), 
r = 94 dB (red, maximum achieved with our experiment), 
r = 100 dB (blue, theoretical limit with our experiment). 


ing the whole sensor, including the sensor head and the 
passive isolation platform, thanks to an excitation plat¬ 
form actuated by PZT. A first test was made by tuning 
the minus-K vertical resonance at a frequency of 2.08 
Hz, identical to the excitation frequency. This setup al¬ 
lows to obtain a pure sinusoidal excitation at a given 
frequency. The 2.08 Hz excitation frequency is chosen to 
stand within the interferometer bandwidth (0 —>■ 1/2T ~ 
0 —>■ 10 Hz), to ensure that the demonstrated rejection 
comes from the differential measurement and not from 
the natural interferometer filter. Moreover this frequency 
is chosen to avoid aliasing effects with the experimental 
repetition rate of 4 Hz. Then, for different vibration am¬ 
plitudes Avib (from few 10“^(; to few 10“^g) the differen¬ 
tial acceleration is derived by fitting ellipses containing 
N = 2000 points. Figure [5] shows three of these ellipses 
for vibration amplitudes of 3.2 x 7.6 x lO^^g and 

2.4 X 10~^g. Between each ellipse acquisition, the fringes 
amplitudes and the population offsets are estimated by 
fitting the interferometric fringes obtained thanks to the 
correlation between the atom and the mechanical ac¬ 
celerometer. 

The results concerning the characterization of the vi¬ 
bration noise rejection are shown in Fig. [6] When the 
amplitude of vibration is lower than typically few rag, the 
atom interferometer sensitivity is still limited by the de¬ 
tection noise and remains nearly constant over this range. 


Conversely, when the amplitude is higher than few mg, 
the vibrations limit the sensitivity and therefore the res¬ 
olution becomes proportional to their amplitude. Under 
these experimental conditions the vibration rejection ra¬ 
tio r, defined as 

r = 20.1og( f-"- ) (14) 

\VN X aAgJ 

can be experimentally evaluated at 88 dB, where cfaq 
is the resolution on the differential acceleration. In ad¬ 
dition to the resolution losses, the vibrations induce a 
drop of contrast through the Doppler shift because of 
additional vertical acceleration. Some contrast losses 
and offset population modifications come from rotations 
(through Coriolis effect [s^, angular and centrifugal ac¬ 
celeration). These rotations are due in our case to a 
non purely vertical excitation with our platform. These 
effects are measured, during the measurement session 
thanks to the mechanical accelerometer and additional 
rotation sensors, and are post-corrected during the dif¬ 
ferential phase extraction from the ellipse {cf. FIGl9]). 
Taking these inertial effects into account, the rejection of 
vibration then reaches a level of 94 dB (a factor 50 000) 
between both isotopes. This vibration rejection level, 
obtained with two different species, is extremely encour¬ 
aging, and chiefly demonstrates the robustness and relia¬ 
bility of a simultaneous differential acceleration measure¬ 
ment with ®^Rb and ®®Rb for testing the WEP in various 
environmental conditions. 

In the next sections we will study what are the limits 
that can explain this level of vibration rejection. 

B. Limit - Wave-vector mismatch 

The first limit comes from the wave-vector mismatch 
between both species. This difference makes both scale 
factors slightly different and any spurious acceleration 
signals are not perfectly canceled by the differential mea¬ 
surement. This is illustrated by terms 7, 8 and 10 in Ta¬ 
ble HI In our experiment 6k/keg = 5 x 10“®, limiting the 
rejection ratio at the same level i.e. 106 dB. This mis¬ 
match is dictated by the way in which the laser lines are 
generated. Indeed, by looking back at the Raman laser 
spectrum (Fig. [S]), the fact that the carrier frequency 
is common to both Raman pairs makes impossible an 
equalization of wave-vectors. 

Figure [7] displays the acceleration transfer function of a 
differential dual-species atom interferometer with a wave- 
vector mismatch (differential 2 curve) compare to the one 
of a single species interferometer (single species curve). 

C. Limit - Difference in temporal response 
functions (Rabi frequencies mismatch) 

Another source of difference between both scale factors 
is due to the difference of temporal response functions, 
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FIG. 7. Normalized acceleration transfer function of a sin¬ 
gle species atom interferometer (black line, single species), a 
dual-species atom interferometer with wave-vector and Rabi 
frequencies mismatches (differential 1, red line), a dual-species 
atom interferometer with a wave-vector mismatch only (dif¬ 
ferential 2, blue line). The wave vector mismatch is Sk/kea = 
5 X 10~®, The Rabi frequencies mismatch between ®^Rb and 
®®Rb is 0.1 (10%). 


^f{t) = fsbit) — /87(0: to accelerations. In our experi¬ 
ment, the atomic beam splitters and mirror are realized 
by the same laser beam for both isotopes. That is why 
the pulse duration, r, and the free evolution time, T, will 
be inherently perfectly matched. Therefore, the response 
function difference comes only from the light-atom inter¬ 
action making the scale factor dependent on the Rabi 
frequency associated to each Raman transition. That is 
why a Rabi frequency mismatch will limit the vibration 
rejection, it is illustrated by term 9 and 10 in Table HI 

In order to emphasize this limitation, we experimen¬ 
tally compared the difference in response function, in the 
low frequency limit, between both isotopes. For this mea¬ 
surement, a very large number of fringes is recorded (~ 
500) for both isotopes by sweeping the frequency chirp a. 
The scale factors are extracted thanks to a sinus fitting, 
and the difference / fsrit) dt — J fsbit) dt = J Sf{t) dt = 
"^eq,s 7 ~ '^eq ,85 accessible. ^ corresponding to 

the temporal part of the scale factor {cf. Eq. (IA5I) in 
appendix lAl) associated to the isotopes i. Figure |S] shows 
the behavior of this difference as a function of the Rabi 
frequency mismatch during the first and third pulses. 
The Rabi frequency mismatch is tuned by changing the 
micro-wave power at 3.035 GHz injected into the phase 
modulator, modifying both Raman pairs powers. The ex¬ 
perimental value of the difference in response functions 
follows the same trend as the theoretical prediction de¬ 
rived from Eq. ([5]) which confirms the validity of the 
predicted scale factor. 

Figure [3 displays the acceleration transfer function of a 


2(Qf-Q“)/(Q*+Qf) 
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FIG. 8. Difference in response functions J fs7{t)dt — 
J fs5{t)dt = f 5 fs 7 {t)dt = Tg,j 87 - Teg,85. m the low fre¬ 
quency limit, between ®^Rb and ®®Rb as a function of the 
Rabi frequency mismatch during the first and third pulses. 
The experimental measurement (black square) is compared 
to the theoretical value (red line) as expressed in Eq. ([5|). 
The horizontal uncertainty corresponds to the uncertainty on 
the determination of the Rabi frequency. The vertical un¬ 
certainty corresponds to the 1-a uncertainty coming from the 
determination of the scale factor by fitting the interferometric 
fringes. 

differential dual-species atom interferometer with a wave- 
vector mismatch to which is added a Rabi frequency mis¬ 
match corresponding to our experimental conditions (dif¬ 
ferential 1 curve). It is very interesting to notice that in 
addition to the rejection loss at low frequency, the in¬ 
terferometer cut off becomes also less efficient at higher 
frequency. In our experiment, the Rabi frequency mis¬ 
match mainly comes from the impact of additional laser 
lines as it has already been explained (c/. section IVl and 
appendix iBl) . In these conditions the vibration rejection 
ratio limitation is evaluated at 100 dB. 


D. Limit - Drop of contrast because of inertial 
effects 

During the vibration rejection experiment, additional 
accelerations and rotations emerge. These inertial effects 
will induce perturbations of the population offsets and a 
drop of the interferometer contrast. They can be visu¬ 
alized in Fig. [5] where the larger the vibrations are the 
noisier the ellipse is. It results a drop of the SNR and thus 
a deterioration of the resolution of the sensor. Neverthe¬ 
less, if the differential atom accelerometer is hybridized 
with others inertial sensors, these spurious inertial effects 
can be measured and post-corrected. The Figure IH shows 
the increase of the SNR after this post-correction. This 
is what has been done to obtain the vibration rejection 
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Normalized Ellipse (vibrations at 2.4x10 g) 
Raw data 
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Normalized Ellipse (vibrations at 2.4x10'^g) 
Post-corrected data 
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FIG. 9. Normalized ellipses obtained with a vibration level 
of 2.4 X 10~^p. At the top, the ellipse is plotted with the raw 
data. At the bottom, the ellipse with a better SNR resulting 
from the post-correction of spurious additional inertial effects. 
This correction is achievable thanks to the hybridization of 
the atom interferometer with classical inertial sensors. 
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ratio of 94 dB (cf. Fig. [51). 

At vibration amplitudes greater than ^ 1 mg, the ac¬ 
celerations in the direction of fceff are large enough that 
the Raman pulses are Doppler shifted out of the Raman 
resonance condition. This leads to an estimated drop of 
contrast of 20 % considering vibrations level of ^ 40 mg 
(peak to peak ; in the following, all the numbers will be 
given according to this level of vibrations). 

Then, these vertical accelerations displace the atoms 
from the ideal detection position by a distance up to 1.2 
mm. The detection efficiency is thus modified in a differ¬ 
ent way for each detection pulse. The normalization of 
atomic populations is also modified which creates fluctu¬ 
ations of the ellipse center up to 5 %. 

The interferometric phase depends on the transverse 
motion of atoms [s^. The rotations can thus induce a 
drop of contrast because of the averaging over the spa¬ 


FIG. 10. Drop of contrast as a function of the amplitude of 
vibrations for ®^Rb (black squares) and for ®®Rb (red circles). 
The drop being induced by spurious inertial effects. 


tial and the velocity spread distributions of the atomic 
clouds. The rotations induce a drop of contrast through 
Coriolis effect [l^, estimated at almost 90 %, for a tem¬ 
perature of the atomic clouds of 2 /rK, which is the most 
disruptive inertial effect. Moreover, angular accelera¬ 
tions, respectively centrifugal accelerations, may create 
additional phase shifts, as a function of the position of the 
atoms, which are averaged over the spatial distribution 
(of about 3 mm FWHM). The resulting drop of contrast 
is estimated at about 60 %, respectively negligible. 

Figure flOl presents the experimental drop of contrast 
induced by all spurious inertial effects previously men¬ 
tioned for both isotopes of rubidium. At low vibration 
amplitude, the ®^Rb contrast is lower than the ®^Rb con¬ 
trast. This is explained by the larger impact of additional 
laser lines on ®^Rb and by its higher temperature. The 
same relation between the contrast is also observed at 
higher vibration level. This is also mainly explained by 
the higher temperature of ®®Rb making the drop of con¬ 
trast due to the Coriolis effect and angular accelerations 
stronger than for ®^Rb. 

In conclusion on the vibration noise rejection and its 
limitations, the rejection ratio of 94 dB (corresponding 
to a factor 50 000) experimentally measured is close to 
the theoretical limits estimated at 100 dB. This demon¬ 
strates that our differential atom interferometer works 
very close to its ultimate performances concerning the 
common-mode vibration noise rejection. The remaining 
difference is attributed to the non-perfect correction of 
the additional inertial effects induced by the excitation. 
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VIII. CONCLUSION 


This experiment is intended to experimentally study 
some limits which could affect a quantum based test of 
the WEP by atom interferometry. 

In this paper we have focused on the sensitivity and 
the resolution of our dual-species atom interferometer for 
measuring the differential acceleration and the profits on 
the common-mode vibration noise rejection. We did not 
study the accuracy of the measurement by estimating the 
systematic effects. Nevertheless, the theoretical study 
about the differential phase (c/. Table |T| allows to state 
about the intrinsic phase shift inherent to the differential 
measurement by atom interferometry. 

The free-fall accelerations of both test bodies are si¬ 
multaneously compared and leads to a WEP test by ex¬ 
tracting the Eotvos ratio 


_ n 087 - 085 „ 

-T = 2- 


Aa 


(as7 + Oss) (o87 + 085) 


(15) 


The WEP is violated when 7787 ^ 1 , ^ 0 that is why 

the differential acceleration Aa needs to be measured in 
an extremely accurate and sensitive way. 

Concerning the accuracy, the WEP violation signal, 
mainly the term 1 in Table U is aimed to be detected. 
All additional phase shifts larger than this signal, com¬ 
ing from the mutual constant acceleration a because of 
different scale factors between both isotopes (terms 5 to 
7), from the atom velocity (terms 11 to 14), from ex¬ 
perimental systematics effects (term 15, as light shifts, 
Zeeman effects, wave-front aberrations, gravity gradient, 
self-gravity, magnetic fields, mean field effects ...) must 
be accurately known and/or estimated to be canceled. 

Concerning atom velocities, when Rabi frequencies are 
not equal between the first and the last light pulses (for 
instance because of the cloud expansion in a given lon¬ 
gitudinal laser intensity profile) the dual-species interfer¬ 
ometer is sensitive to the velocity difference between both 
species (term 11 to 14 and mainly term 11). By nam¬ 
ing Sfl the Rabi frequency mismatch between the light 
pulses, the condition Term 11 < Term 1 leads to the fol¬ 
lowing condition on the velocity difference between both 
isotopes : Av < ■ By considering Aa = 10“^"^ 


^-2 


m.s^", T = 1 s, r = 100 /iS, = 10“^ we obtain 

Av < 2 X 10“® m.s“^. This condition is experimentally 
challenging but remains below the conditions on the dif¬ 
ferential velocity targeted by STE-QUEST to counteract 
effects as the gravity gradient or the Coriolis acceleration 


Concerning the sensitivity, future experiments dedi¬ 
cated to test the WEP aim to be limited by the quantum 
projection noise during the atomic populations detection 
|5^ [ 5 ^ . Such working regimes are currently reached 
in state of the art quantum inertial sensors. The main 
source of noise, which is inherently linked to the accel¬ 
eration measurement, and which could thus impair the 
sensitivity is obviously the vibration noise, i.e. a terms 


in TableHl As it has been previously explained, a simulta¬ 
neous handling of both quantum test bodies is required 
to benefit from a common-mode vibration noise rejec¬ 
tion. The use of both isotopes of the same species is 
here a strong advantage because it allows to take advan¬ 
tage of high rejection levels (as demonstrated in this pa¬ 
per) which ensure the high reliability of the differential 
accelerometer in a large range of environmental condi¬ 
tions. On ground, this allows to release some techno¬ 
logical constraints on the isolating vibration systems for 
instance. In micro-gravity environments, according to 
terms 8 and 9, the use of ®^Rb and ®®Rb allows to work 
with white noise levels, in the interferometer bandwidth, 
up to 2 X 10“'^ Tn.s~'^ / \/Hz by considering a relative dif¬ 
ference of 10“® on wave vectors and 10“^ on Rabi fre¬ 
quencies, or even up to 1 x 10“^ /\/Hz if a 10“"* 

level on the Rabi frequency match can be achieved (the 
values ot k, T, T has been here chosen according to STE- 
QUEST parameters This is a strong argument 

for the use of these two isotopes for testing the WEP in 
the international space station (ISS) or in a satellite. In¬ 
deed the vibration noise in the ISS is quite high because 
of human activity, whereas in a satellite, the vibration 
levels remain little-known and strongly depend on the 
satellite itself and on all the instruments that are loaded 
in there. The use of ®^Rb and ®®Rb thus offers a strong 
safety in order to counteract a large number of spurious 
and unexpected environmental effects. 

We have reported here the realization and the charac¬ 
terization of a simultaneous dual-species atom interfer¬ 
ometer which handles both ®^Rb and ®^Rb. We care¬ 
fully studied the expression of the differential phase and 
the measurement scale factor, thanks to the sensitivity 
response function formalism, and explained the method 
used to derive the differential phase. After having tackled 
the main limitations of our previous device (the impact of 
additional laser lines and the non-linearities in the detec¬ 
tion scheme) we showed that our differential accelerom¬ 
eter currently reaches a sensitivity of 1.23 x 10~^g at 1 s 
and a resolution of 2 x 10 “®g after an integration time of 
few hours only. These results have been obtained with a 
small and compact atom sensor. The simultaneous dif¬ 
ferential measurement allowed us to exhibit a vibration 
rejection factor of 50 000 which proves the reliability of 
such a sensor in a large range of environmental conditions 
and its utility for future projects aiming to test the WEP 
with matter-waves. A next step to improve the current 
results could be to implement a new detection scheme for 
better counting the atomic populations of both species. 
Finally, it is necessary to work further on the correction 
of the systematic effects related to the measurement in 
order to reach an accuracy on 77 at a level of 10 “®. 
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mental setup. 


The interferometric phase A<i> can then be derived by 
the formula 


Appendix A: Response function of a Mach-Zehnder 
type atom interferometer to acceleration 

Initially developed for atomic clocks [s^, the sensitiv¬ 
ity function formalism is well adapted to evaluate the 
response of the interferometer to laser phase fluctuations 
and thus to acceleration. 

Considering an infinitesimal phase step 54> of the Ra¬ 
man laser phase occurring at time t the sensitivity func¬ 
tion is defined as in Eq.(l5|). We make the same assump¬ 
tions than in (45l| . i.e. the laser waves are considered as 
pure plane waves, the Rabi frequencies are constant dur¬ 
ing a pulse (square pulses) and the resonance condition is 
fulfilled. Moreover we have extended this framework by 
assuming that the Rabi frequencies are not necessarily 
equal between each pulse (fli ^ ^2 ^ ^ 3 )- We followed 
the same matrix calculation approach than in . With 
a time origin in the middle of the interferometer, the ex¬ 
pression of the sensitivity function is given by 


A$ 



9s{t) 


dt 


dt 


(A2) 


9 sit) 


' 0 , 

sin[ni(t+T+2r)] 

sin(Qir) 

-1, 

sin(n20 

sin(f22T) ’ 

1, 

sin[Q 3 (j:—T—2r)] 
sin(Q3T) 


t < -T- 2r 
-T -2t <t<-T 

— T — T < t < —T 

— T <t<T 

T < t < T + T 
T + T<t<T + 2T 


t>T + 2t 

(Al) 


J 


where ip{t) is the phase of the Raman laser seen by the 
atoms. In the frame associated to the free falling atom, 
this phase is equal to (pit) = keS-r{t), with fit) being the 
position of the Raman mirror compared to the atom. The 
sensitivity function corresponds to the response of the in¬ 
terferometer to the velocity of the free falling atom. In¬ 
tegrating by parts the Ea. (IA2() . it can be shown that the 
interferometric phase can be expressed as in Eq.Q. The 
response function of the interferometer to acceleration is 
thus the primitive integral of the sensitivity function (c/. 
Eq.©) whose complete expression is 


_ tan {^) - ^ tan , 

Oi sin(nir) (cos(nir) - cos + T + 2t)]) + jA tan , 

t + T + T + -^ tan (^) , 

fit) = a, sinWr) (cos(n2t) - cos(fl2T)) + T + ^ tan (%i) , 

-t + T + T+ ^tan (^) , 

a 3 sin( 03 r) (1 - cos W - T - 2t)]) , 

lo. 


t<-T-2T 
-T -2 t <t<-T -T 

— T — T < t < —T 

— T <t<T 

T < t < T + T 
T + T<t<T + 2T 
t>T + 2T 


(A3) 


This triangle-shaped function (c/. FIGH]) evolves as 
cos(Oit) during the pulse i and is maximum in the mid¬ 
dle of the interferometer where the separation between 
both arms of the interferometer is the largest. The par¬ 
ity break of the interferometer symmetry between the 
first and the third Raman pulse makes it sensitive to the 


atom velocity. That is why f{t) is not perfectly zero for 
t<-T- 2t. 

Thanks to this formalism, the expression of the scale 
factor S impacting the measurement of a constant accel¬ 
eration can be estimated by taken into account the finite 
duration of the Raman pulses : 
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s = keS gs (t) t dt 


where E'down and E^p corresponds to the electric field 
respectively going downward and upward and 4> to the 
depth of modulation. 


^eff 

(T + 2t) (t + ^ tan (^) + tan 

+ (?|-n?-nf) 


-T 2 


0 


cotan(r22T) cotan(r2iT) cotan(r23r) 


^2 


Qi 


Qs 


(A4) 

The temporal part of the scale factor is thus given 
at first order by : 


The next step is to describe the probability amplitude 
to go from state |a) to state \b) when this electric field 
is turned on. We consider only a two-photon Raman 
transition with counter-propagating beams (the others 
co-propagating and opposite counter-propagating being 
neglected because of the Doppler shift induced by the free 
fall of atoms). Each component E{n) of the electric field 
couples the two states with an associated Rabi frequency 
rin proportional to 


^eq ■ 






-®down(^ + l)^up(fT-) 


(B2) 


A — nuJG 

where A = ujl— ojg- The resulting “overall” two-photon 
Rabi frequency is thus given by 


Appendix B: Impact of additional laser lines 
generated by modulation on the atom interferometer 

See also reference for a more comprehensive ap¬ 
proach. 

Let us consider an atom corresponding to a A-type 
three-level system with two ground states |a) and |6) sep¬ 
arated by an energy Tiuig and an excited state |e) sepa¬ 
rated by hcoo from state |6). The atom interacts with a 
laser which is retro-reflected by a mirror at a position zm 
from the atom. The laser spectrum is composed of laser 
lines separated in frequency by log and centered around 
the laser carrier pulsation lol. This kind of spectrum 
is obtained with a phase modulated laser at frequency 
LOG (c/. spectrum corresponding to one species in Figl3]). 
The electric field experienced by the atom can be decom¬ 
posed on the basis thanks to the Bessel functions 

Jn, and written as the sum of 

L^down = EoE^^^ i:l=-oo + c.c 

E^p = + c.c 

(Bl) 


+ 00 

^|a)^|b) — ^ ^ 

n— — oo 


(B3) 


In this formalism we do not consider the external state 
of the atom, the degeneracy of the quantum states due 
to the slight difference in momentum h/S.k = 2‘hLOGlc. 
between the different Raman pairs is neglected. This 
approximation is a bit rough but allows to simply explain 
the problem.The overall Rabi frequency is thus expressed 
as 


+ 00 

f^!a)->-|b) OC ^ 

n= —oo 


^inAkzM 

A — ULOG 


+ c.c 


(B4) 


The terms in this sum interfere, which has the effect of 
spatially modulate the Rabi frequency. The global inter¬ 
ferometer contrast will depend on the distance between 
the atoms and the mirror (zm) and thus also on the time 
when each Raman pulse occurs. 
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